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Abstract 

We study the multifractal properties of the current distribution of the three- 
dimensional random resistor network at the percolation threshold. For lattices 
ranging m size from 8^ to 80^ we measure the second, fourth and sixth mo- 
ments of the current distribution, finding e.g. that t/u = 2.282(5) where t is 

the conductivity exponent and v is the correlation length exponent. 
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I. INTRODUCTION 



Quite surprisingly, it was found in the mid-eighties that dynamic phenomena on fractal 
structures often were controlled not by one or two relevant length scales, but rather an 
infinite hierarchy of such length scales [Q. One of the prime examples of this phenomenon 
is the current distribution in the random resistor network at the percolation threshold . 
In spite of the large effort which was invested to understand how such an infinite hierarchy 
could appear, e.g. through studying hierarchical structures yielding to analytic calculations 
no satisfactory general explanation was found. The large numerical effort that was 
invested at that time on the random resistor network focussed on two dimensions. The 
reason for this was that three-dimensional networks were essentially out of reach of the 
computational power available at the time. The aim of the present work is to establish 
that, as in two dimensions, the current distribution in three dimensions is multifractal and 
to determine the corresponding scaling exponents with high presicion. In addition to their 
theoretical interest, these results are of importance in making contact with experimental 
studies on three-dimensional conductor- insulator mixtures and microemulsions p|-[13|. 

To achieve this goal, we made use of the latest developments in iterative solvers and massively 
parallel computers. 

In section ^ we present the model and the method of solution. In section |T| we discuss 
the current distribution through the behaviour of the moments and their exponents, in 
addition to examining the statistical fluctuations of the moments. 

II. MODEL AND NUMERICAL SOLUTION 

We consider a three-dimensional cubic lattice of size with periodic boundary condi- 
tions in the x and y directions. For the z direction the boundary conditions are as follows: 
At z = 1 and z = L we place two plates with a constant potential difference set at a value 
of 1. Therefore, the length of the lattice in the z direction is {L — 1), while in the x and 
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y directions the length is L. This geometry was chosen because the data layout becomes 
optimal on the Connection Machine CMS which we used for our computations. 

All bonds are visited and, with a probability p, a resistor is placed. All resistors have the 
same resistance, 1. We set p equal to the bond percolation threshold for the cubic lattice, 
Pc = 0.2488, as determined by Stauffer et al. After all the bonds have been visited, 

a (parallel) cluster finding algorithm is applied to determine if there is a spanning cluster 
that connects the two plates. If no such cluster exists, the procedure is repeated until one 
is obtained. This cluster finding algorithm is very fast but only determines if there is such 
a cluster, not its exact geometry. 

The equations to be solved are the usual current equations (Kirchhoff 's equations) on the 
lattice which can be easily solved using the conjugate gradient or related iterative algorithms 



| T5| . We wrote the program in CM Fortran and used the iterative solvers in the Connection 
Machine Scientific Software Library (CMSSL) which contains thirteen of them. Which one 
to use, depends on the particular matrices one is dealing with. For this problem we found 
that the quasi-minimized CCS (QCGS) [|TB[ algorithm performs very well. 

Our stopping condition is for the residual to be less than 10"^^, which for the biggest 
lattices gave a true accuracy of about 10~^. We estimated this by calculating the conductivity 
of each realization in two ways: 1) Calculate the total current crossing an XY plane at 
z = L/2, which, knowing the potential difference (=1), gives the conductivity, and 2) 
calculate the second moment of the currents. This is again equal the conductivity as the 
externally applied potential difference is 1. On the biggest lattices, for which the accuracy 
at which we can determine the currents is the lowest, they agreed to within 0(10^^). As 
it is too time consuming compared to solving the Kirchhoff equations, we made no effort 
to identify the bonds that did not belong to the spanning cluster. We simply solved for 
the currents keeping the disconnected bonds and dangling ends in the system. This made 
it impossible to get an accurate count of the number of current-carrying bonds. We will, 
therefore, give results for the second, fourth, and sixth momemts of the currents, leaving 
out the zero'th moment. The number of realizations for L = 8, 16, 32, 48, 64, and 80 were 



12000, 9198, 2765, 1120, 1913, and 953 respectively. 



III. CURRENT DISTRIBUTION 

Solving for the currents in the bonds allows us to calculate the nth moment of the current 
distribution, given by 

{nv = ^T.T.^m. (1) 

R k(R) 

where R denotes the realization, Nr is the number of these realizations, and ik(R) is the 
current in the fcth bond in realization R. These moments are calculated in the constant 
voltage ensemble where the potential difference is kept constant realization to realization. 
We expect the moments of the currents to scale as 

{^v ~ Lf-\ (2) 

where = L — 1 is the length of the system in the z direction, and x{n) is the exponent 
of the nth moment. Figures |l|a, ^ and |^ show, on a log- log scale, the second, fourth and 
sixth moments respectively. The second moment measures the total dissipated power, and 
since the externally applied voltage difference is 1, it directly gives the conductance of the 
network. So, Figure [l|a gives the scaling of the conductance as a function of L^, giving the 
exponent x{2) = 1.282(5). Figure [pD shows the scaling of the variance of the conductance 
distribution as a function of L^. We see that it scales with the same exponent as the 
conductance itself. The same is true for the other moments. Thus, the relative fluctuations 
neither grow nor decrease with lattice size. The scaling of the conductivity with is 
obtained by simply dividing the results for the second moment by Lz (since this is the 
conductance) giving t/v = 2.282(5), where t is the conductivity exponent and v is the 
corrrelation function exponent. Our value for t/v is in agreement with, but more precise 
than, the value determined by Gingold and Lobb ||T^ who obtained tjv = 2.276(12). With 
V = 0.88, we therefore have that t = 2.01. The experimental values reported for this 
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exponent ranges from 1.2 to 2.1 for the measurements based on microemulsions P-[T3 
and for the measurements based on conductor-insulator mixtures, the values 2.0 ± 0.2 0, 
1.85 ±0.25 0], and 1.6 ±0.1 §. 

The exponent of the fourth moment is related to the scaling of the Nyquist noise of the 
random resistor network through the fluctuation-dissipation theorem as shown by Rammal 
et al. [Q. Thus, as for the second moment, the fourth moment is related to a macroscopic 
quantity and therefore is of direct experimental interest. However, no relation between 
the sixth moment of the current distribution and a directly measurable quantity has been 
identified. We determine a;(4) = 3.920(6) = and x(6) = 6.477(10). 

If y{n) is the exponent of the nth moment in the constant current ensemble, we have 
(^")c ~ -L^*-"-*. It is then easy to show that y{n) = x{n) — nx(2). This way, we can easily 
change between the two ensembles. Figure ^ shows y{n) as a function of n, where we see 
that it behaves in the classic multifractal way PJT8[|: The three exponents we show do not 
fall on a straight line and as n increases they approach a constant that must equal l/u P,p!9|. 
This is shown by the dashed line in the figure. To show this relation even more clearly, we 
plot in Figure |^ y{n) versus 1/n, and where we use for y{oo) the best value we found for 
u = 0.88, and references therein. The dashed line is merely to guide the eye. It connects 
the point at 1/n = to that at 1/n = 0.5. We see that the exponents for the fourth and 
sixth moments are in agreement with this line although the sixth moment is starting to lose 
precision. These y{n) plots demonstrate the multifractality of the current distribution in 
this network, but the convergence to 1/z/ as n ^ oo is rather slow, especially when compared 



to the two-dimensional case 18 



The sample to sample fluctuations of the values of the conductivity (G) and the fourth 
moment (related the Nyquist noise strength, let us call it i^' = {i^)v) also yield interesting 
information about the system ||21|. We have demonstrated in Figure |I]b that the variance 
of the conductance (from sample-to-sample fluctuations) scales as the conductance itself 
with respect to L^. Exactly the same behaviour is observed for the variances of the higher 
moments. Therefore, the distributions of G/ (G) from different size systems will collapse onto 



a single distribution. The same is true for K/{K). To characterize these distributions, we 
examined on semi-log scale these distributions against {G—{G))^/ {G)^ and {K—{K))^/ {K)"^, 
respectively. Such plots should yield straight lines for Gaussian distributions. This way we 
found the distributions not to be normal, and a similar procedure showed them not to 
be lognormal. In Figures ^ and |^, we show the distributions of G/{G) and K/{K) for 
different lattice sizes on semi-log scale plotted against G/{G) and K/{K). For the larger 
values, we find in both plots straight lines which indicate exponetial distributions of the 
form N{G/{G)) ~ exp(-2G/(G)) and N{K/{K)) ~ exp{-1.3K/{K)), respectively. 



IV. CONCLUSIONS 

By using a combination of efficient algorithms and a massively parallel computer, we 
were able to do a high precision study of the distribution and moments of currents in a three 
dimensional network at the percolation threshold. 

Our result for the conductivity exponent is in agreement with but more precise than 
previous values. In addition we evaluated the exponents for the fourth (related to 1/f and 
Nyquist noise) and sixth moments. The values we have found support the notion of a 
multifractal current distribution for the three dimensional network. 

We have furthermore studied the sample-to-sample fluctuations and the distributions 
of the second (conductivity) and fourth (noise) moment of the current distribution. We 
find that the relative fiuctuations scale as the moments themselves with lattice size, and 
that the underlying statistical distributions appear to be exponential rather than Gaussian. 
Furthermore, the distribution of conductivities at the percolation threshold has proven to 
be an important ingredient in the formulation of scaling theories for the optical properties 
(ac conductivity, refiectivity, transmittivity) of two-dimensional systems such as semicon- 



ductor metal films |^ and metal-insulator composites ||23[ . Similar scaling theories in three 
dimensions, when constructed, will also need the distribution of the conductivities for three- 
dimensional systems at the percolation threshold, which we have presented in this paper. 
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FIGURES 

FIG. 1. a) The second moment of the current {i.e. the conductance G) scales as a power 
of the system size Lz, with an exponent equal to —1.282(5). b) The variance of the conductance 
distribution scales exactly as the conductance itself with the system size, Lz- 

FIG. 2. Power law scaling of the fourth moment as a function of the system size Lz. The 
exponent is —3.920(6). 

FIG. 3. Power law scaling of the sixth moment as a function of the system size L^. The 
exponent is -6.477(10). 

FIG. 4. y{n), the exponent of the nth moment in the constant current ensemble, as a function 
of n. As n increases, y{n) approaches the dashed line which is given by l/v. 

FIG. 5. Same as Figure 4, but shows y{n) as a function of 1/n. This shows very clearly that 
as n — > DO, y(n) ^ l/i/. The dashed line is a guide to the eye and simply connects the points at 
1/n = and 1/n = 0.5. 

FIG. 6. Histogram of G/ (G) on semilog scale. It shows that the distribution for larger values 
of G/{G) is exponential. 

FIG. 7. Histogram of K/ (K) on semilog scale. It shows that the distribution for larger values 
of K/ [K) is exponential. 
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Figure 4 

G.G. Batrouni, A. Hansen and B. Larson 
' 'Current Distribution in the 3D Random Resistor Network" 



1.35 



1.30 - 



1.25 - 



c 



1.20 - 



1.15 - 



1.10 



0.00 




Figure 5 
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'Current Distribution in the 3D Random Resistor Network" 
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